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Abstract 

Using the first-principles GW-Bethe-Salpeter equation method, here we study the excited-state 
properties, including quasi-particle band structures and optical spectra, of phosphorene, a two- 
dimensional (2D) atomic layer of black phosphorus. The quasi-particle band gap of monolayer 
phosphorene is 2.15 eV and its optical gap is 1.6 eV, which is suitable for excitonic thin film solar 
cell applications. Next, this potential application is analysed by considering type-II heterostruc¬ 
tures with single layered phosphorene and transition metal dichalcogenides (TMDs). These het¬ 
erojunctions have a potential maximum power conversion efficiency of up to 12%, which can be 
further enhanced to 20% by strain engineering. Our results show that phosphorene is not only a 
promising new material for use in nanoscale electronics, but also in optoelectronics. 
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I. INTRODUCTION 


Q 


The discovery of graphene in 2004 [jj was a significant breakthrough in materials science 
and since then there has been sustained research interest in graphene and a rush to discover 
other stable two-dimensional (2D) materials. 2D materials have the thickness of one or a few 
atomic layers and have markedly different material properties than their bulk counterparts 
due to the quantum conhnement effect. Since graphene, there have been several advances 
in the field of 2D materials such as the discovery of semiconducting monolayer transition 
metal dichalcogenides (TMDs) { 2 , 3|. The growing number and variety of 2D materials has 
fuelled interest in the use of 2D materials in the development of novel nanoscale devices. 

A recent development is the experimental isolation of a single layer of bulk black phos¬ 
phorus, also known as pho^horene 4l-l23l|. Like graphene, phosphorene was first obtained 
by mechanical exfoliation |7|, l8|. Liquid exfoliation, which is a scalable process, has also 
been demonstrated as a possible alternative means of producing phosphorene jd]. Several 
experiments have demonstrated that phosphorene is a direct-gap semiconductor and also 
has a high hole mobilityj?, 24|. These characteristics make phosphorene attractive for 
use in electronic and optoelectronic devices 2^. Field effect transistors (FET) based on 
few-layer phosphorene were shown to have high on/off ratios [^. It also has the potential 
to be used as an anode material in lithium-ion batteries . Another potential application 
is in thin him excitonic solar cells as phosphorene has a predicted band gap in the visi¬ 
ble region 26|]. Excitonic solar cells (XSCs) based on some 2D materials, such as M 0 S 2 , 


WSe 2 , graphene, h-BN, SiC 2 anc 
eration of thin him solar cells 


anc 

bil 

26 

31 


ayer phosphorene, are potentially seen as the new gen- 


31| . and they might have higher efficiencies than existing 


XSCs which typically have less than 10% efficiency 3^. Till now, despite the limitations 
of fabrication methodologies for such 2D solar cells, there has been some progress on the 
fabrication of few-layer heterostructures such as graphene-WS2[33|, [3^, graphene-MoSo [35 1 
and phosphorene-MoS 2 |. 

In this paper we study the excited-state properties of monolayer phosphorene, and eval¬ 
uate the viability of monolayer phosphorene as one building block of of an excitonic solar 
cell heterostructure. For the other building-block material in the heterostructure, the semi¬ 
conducting monolayer TMDs, which have been extensively researched, are chosen. The 
semiconducting TMDs considered here include semiconducting M 0 S 2 , MoSe 2 , MoTe 2 , WS 2 , 
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FIG. 1. (Color online) GW-calculated (solid lines) and HSEOG-calculated (dash lines) band struc¬ 
ture of monolayer phosphorene. The Fermi level is set on the top of the valance band. The ball and 
stick model of phosphorene (light color of the upper layer) with primitive cell and the reciprocal 
lattice with the high symmetry points are shown in inset. 

WSe 2 , WTe 2 , TiS 2 and ZrS 2 . This paper is organized as follows: In Sec. II, we first intro¬ 
duce the structures of geometrically optimized monolayer phosphorene and computational 
details. Sec. Ill is the main part of results and discussions, which includes three sub¬ 
sections: In SubSec. A, we show the excited-state properties of monolayer phosphorene, 
including in quasi-particle band structures and optical spectra. In SubSec. B, we calculate 
the excited-state properties of 8 semiconducting TMDs. In SubSec. C, the band alignment 
of phosphorene and TMDs is presented. Then the power conversion efficiencies (PCE) of 
the phosphorene-TMD heterostructures are discussed. In Sec. IV, we conclude our studies. 


II. STRUCTURES AND COMPUTATIONAL DETAILS 

Phosphorene has a puckered honeycomb structure as shown in inset of Figure 1. The 
underlying lattice is rectangular which leads to anisotropy in the band structure and optical 
properties. The calculated lattice constants are a = 4.58 A and b = 3.30 A which are consis- 
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tent with results in literature SL The calculations were performed using the Vienna Ab-initio 
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401 1 . and a projector augmented wave (PAW) basis set was 


42| . Geometry optimisation was done using density function theory (DPT) with 


the generalised gradient approximation (GGA) using the Perdew-Burke-Ernzerhof (PBE) 
exchange correlation functional. A 14 x 10 x 1 /c-point grid was used for phosphorene while 
anllxllxl /c-point grid was used for TMDs. The geometry was relaxed until the force 
acting on the atoms was less than 0.01 eV/atom. To ensure that the interlayer interaction 
is negligible, the out of plane lattice parameter, which is perpendicular to the plane of the 
material, was set as at least 15 A. Band gaps and band structures were calculated using the 
GW method. The band gaps/structures were also calculated through the screened exchange 
hybrid density functional by Heyd, Scuseria, and Ernzerhof (HSE06) for reference. A high 
number of empty conductions bands is necessary for the convergence of the absolute band 
positions!^. A total of 1024 bands were used for phosphorene and 1536 bands were used 


for TMDs. Single shot GW calculations, GqWo, were performed on phosphorene to obtain 
the ground state band energies. For TMDs, one eigenvalue u pda te is performed to obtain 
the expected direct band gap of trigonal prismatic TMDs 44l-l47l|. The band structure was 
then interpolated from Wannier functions rather than evaluated directly at discrete /c-points. 
This was done using the WANNIER90 library and the VASP2WANNIER interface. The 
optical gap was calculated by solving the Bethe-Salpeter equation (BSE). The GW method 
has more computational demands and the band gaps converge with a much smaller number 
of bands 431]. Thus to streamline optical calculations, the GqWq calculation was run with 
128 bands for phosphorene and 192 bands for TMDs before solving the BSE. This produced 
the frequency dependent dielectric tensor that was used to calculate the absorption spec¬ 
trum. By solving the BSE, electron-hole interactions such as excitons are accounted for in 
the dielectric tensor. For phosphorene the x- and y-components of the dielectric tensor were 
treated separately because of the anisotropy but for TMDs, the average value of the x- and 
y-components was used due to crystal symmetry. 
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FIG. 2. (Color online) GW-calculated (solid lines) and HSEOG-calculated (dashed lines) absorption 
spectrum of phosphorene for armchair polarised light. The energy values of first absorption peak, 
observed in the experiment j3], is presented by a black solid circle for reference. 

III. RESULTS AND DISCUSSIONS 


A. Excited-state properties of phosphorene 


The calculated band structure of phosphorene is shown in Figure 1. The GW calculation 
predicts a band gap value of 2.15 eV and the gap is approximately located at the T point, 
which is much higher than the DFT-calculated band gap of 0.91 eV 9|. The conduction 
band minimum (CBM) and the valence band maximum (VBM) are not exactly aligned in 
the GW-calculated band structure, but they are sufficiently close to be considered as a direct 
band gap |^. A similar computational approach W Tran et al. predicted a band gap of 2.0 
eV and a comparable profile of the band structure [l8j. The CBM position is around -4.25 eV 
with respect the the vacuum level. Besides the GW method, the hybrid density functional 
HSE06 was used to calculate the band structures and band gap for reference in Fig. 1. As 
can be seen, the HSE06-calculated band gap is 0.6 eV lower than the GW-calculated gap, 
which is in agreement with previous calculations jo]. 
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FIG. 3. (Color online) (a) The band structures of M 0 S 2 with a ball and stick model of the trigonal 
prismatic geometry, (b) The band structures of TiS 2 with a ball and stick model of the octahedral 
geometry. The reciprocal lattice with high symmetry points is shown in the inset. The Fermi level 
is set on the top of the valance band. 


The optical gap of phosphorene is calculated using the GW-BSE approach and is de¬ 
termined to be 1.6 eV. This is the first optical peak of the absorption spectrum for light 
polarised along the armchair direction as seen in Fig. 2. This optical gap is slightly larrer 
than the experimental photoluminescence measurement value of 1.45 eV of phosphorene |8|, 
but smaller than the HSE06-calculated value of 1.8 eV. The optical gap of 1.6 eV is much 
lower than the electronic band gap of 2.15 eV which suggests that significant excitonic effects 
are present in phosphorene. The exciton binding energy of 0.55 eV in monolayer phospho- 
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TABLE I. Summary of the lattice constants a and b, the band gap Eg, valence band maximum 
Evbm, conduction band minimum Ecbm and optical gap Eopt 


Material 

Lattice 

Structure 

a (A) h (A) E 

^.g (eV) Evbm (eV) Ecbm (eV) E 

opt (eV) 

Phosphorene Rectangular 

- 

4.58 3.30 

2.15 

-6.40 

-4.25 

1.6 

M 0 S 2 

Hexagonal 

Trigonal Prismatic 

3.16 - 

2.68 

-6.57 

-3.89 

2.3 

MoSe2 

Hexagonal 

Trigonal Prismatic 

3.29 - 

2.39 

-5.97 

-3.58 

2.2 

MoTe2 

Hexagonal 

Trigonal Prismatic 

3.52 - 

1.74 

-5.43 

-3.69 

1.7 

WS 2 

Hexagonal 

Trigonal Prismatic 

3.16 - 

2.94 

-6.47 

-3.53 

2.5 

WSe2 

Hexagonal 

Trigonal Prismatic 

3.26 - 

2.70 

-5.84 

-3.14 

2.6 

WTe2 

Hexagonal 

Trigonal Prismatic 

3.52 

1.98 

-5.35 

-3.37 

1.9 

TiS2 

Hexagonal 

Octahedral 

3.37 - 

1.88 

-6.80 

-4.88 

- 

ZrS2 

Hexagonal 

Octahedral 

3.58 - 

2.65 

-7.56 

-4.81 

- 


rene is quite huge. Both the self-energy corrected large electronic gap and small optical gap 
of phosphorene indicate significant many-electron effect in phosphorene. 


B. Excited-state properties of TMDs 

Using a similar approach to the above, the geometry, band positions and band structures 
of 8 semiconducting TMDs are calculated. A summary of the lattice constants, absolute po¬ 
sitions of the valence band maximum (VBM) and conduction band minimum for the different 
materials is shown in Table 1. The TMDs are categorised into either trigonal prismatic or 
octahedral TMDs, both of which have a hexagonal lattice but different coordination of the 
atoms within the unit cell. The two different structures are shown in the inset of Fig. 3. 
For trigonal prismatic TMDs, the lattice constant is largely determined by the size chalcogen 
atom which increases with atomic number (sulphur to selenium to tellurium). Furthermore, 
as can be seen, phosphorene has a large exciton binding energy compared to TMDs because 
of its unique quasi ID band dispersions. 

The trigonal prismatic TMDs have a direct band gap the the K point [see Supplemen¬ 
tary Materials]. They have the similar band structures. The band structure of M 0 S 2 is 
shown in Fig. 3(a) as an example. For trigonal prismatic TMDs, the band gap decreases 
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with increasing chalcogen atomic number. The CBM position also decreases with increas¬ 


ing chalcogen atomic number. These trenc 
other similar studies of these materials 43 


s and the band gap values are consistent with 


44j. kOml]. Octahedral TMDs have an indirect 


band gap between the T and M points [see Supplementary Materials]. This is also in 
agreement with band prohles calculated in previous studies, 
structure of TiS 2 is shown in Figure 3(b). 
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FIG. 4. (Color online) Band alignment of phosphorene with semiconducting monolayer TMDs. 


Note that many important properties and potential device applications of semiconductors 
are not determined entirely by the band gap only. The band alignment and corresponding 
band offsets (the relative band-edge energies) of two or more semiconductors are other 
fundamental/critical parameters in the design of heterojunction devices 


43 


example, the 2D heterostructure devices 


effect transistors 


or photocatalytic water splitting 


50 


43 


m 


54; 56], for 


54l |. field 


49| and p-n diodes 3^. Chemical trends of the band offesets provide 


a useful tool for predicting catalytic ability of TMDs-based heterojunctions. Figure 4 
shows the band alignment (using the vacuum level as reference) of phosphorene with 8 
semiconducting monolayer TMDs. It can be seen that the CBM of trigonal prismatic TMDs 
is higher than that of monolayer phosphorene, while the CBM of octahedral TMDs is lower 
than phosphorene. Thus, trigonal prismatic TMDs could function as the donor whereas 




























































octahedral TMDs would function as the acceptor in heterostructure with phosphorene. The 
optical gap of 8 TMDs is calculated, which is determined from the absorption spectrum. 
In most cases, the optical gap is slightly lower than the electronic band gap as shown in 
Table 1. Notice that because of strong many-electron effects in some 2D materials, such 
as phosphorene and M 0 S 2 , we might not get the accurate band offset parameters in some 
cases without considering the excited-state effect. For example, the HSE-calculated CBM 


19| (-4.25 eV 50|) and phosphorene is -3.94 eV 


band energy of monolayer M 0 S 2 is -4.21 eV 

1^ (-3.92 eV 1^1). Thus, phoshporene (M 0 S 2 ) is the donor (acceptor) in the phosphorene- 
M 0 S 2 2D heterojunction [l9| and the predicted PCE is 17.5% 19|. However, if we consider 
many-electron effects, the GW-calculated CBM band energy of monolayer M 0 S 2 is -3.89 eV 


(-3.74 eV 


43|) and phosphorene is -4.25 eV. Interestingly, phoshporene (M 0 S 2 ) becomes the 


acceptor (donor) instead and the predicted PCE is reducing to 10%. Regarding to other 2D 
materials without strong many-electron effects, both HSE and GW can give similar band 


offset 


43 


50|. 


A model developed by Scharber et ah for organic solar cells (57| and later adapted for 


exciton based 2D solar cells 


26 


29| is used to predict the maximum PCE, rj, based on the 


hll-factor, (3ff, open circuit voltage, Voc, and short circuit current, Jg. 

f^FpVocJsc 


V = 




( 1 ) 


solar 


where Psoiar is the total incident solar power per unit area based on the Air Mass (AM) 1.5 
solar spectrum 58|, l59| . The £11 factor is the ratio of power output at the maximum power 
point to the product of the open circuit voltage and the short circuit current. The fill factor 
is estimated to be 0.65 from literature. The Voc, in units of V, and Jgc, in nnits of A/m^, 
are estimated in the limit of 100% external quantum efficiency as 


Voc — - {Eopt — ^Ecbm — 0 . 3 ) ( 2 ) 

Jsc = e f (3) 

where e is the elementary charge, is the donor optical gap, AEcbm is the conduction 
band offset and P{hu) is AM 1.5 solar spectrum. In equation [2], the constant 0.3 eV is an 
empirical parameter that estimates losses due to energy conversion kinetics. 

Using this model the TMDs are paired with phosphorene. The material with the lower 
CBM is the acceptor and the material with the higher CBM is the donor. Phosphorene 
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FIG. 5. (Color online) Power conversion efficiency of potential excitonic thin film solar cell hetero¬ 
junctions. The labels indicate the material that complements phosphorene. 


is the donor when paired with both octahedral TMDs and the acceptor when paired with 
trigonal prismatic TMDs. The maximum PCE values for these eight heterostructures are 
marked on Fig. 5. Of the eight heterostructures, phosphorene-ZrS 2 and MoTe 2 -phosphorene 
have the highest PCE value of 12%. This efficiency is higher than that achieved by existing 
excitonic solar cells, and comparable to the proposed 2D g-SiC 2 /GaN (14.2%), PCBM/CBN 
(10-20%) Q, and bilayer-phosphorene/MoS 2 (16-18%) solar cells. 


Actually, a further observation based on this model is that a solar cell with a phosphorene 
donor could have maximum PCE values of up to about 20% with an appropriate choice for 
the acceptor. The conduction band offset (CBO) between phosphorene and TiS 2 is 0.63 
eV while the CBO between phosphorene and ZrS 2 is 0.56 eV for the two cases here where 
monolayer phosphorene is the donor. This translates to a signihcant drop in Voc in the 
model which results in a lower PCE. In addition to looking out for new materials that 
have a better band alignment with phosphorene, means of tuning the properties of both 
the donor and acceptor can be considered. For example, strained phosphorene may be 


used as the acceptor materia 
band structure of materials 


because the strain effect is a well-known method to tune the 


15 


, ll6|]. The band structures of different strained monolayer 
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FIG. 6. (Color online) (a) The band structure of phosphorene under applied 2% armchair com¬ 
pressed strain. The Fermi level is set on the top of the valance band. The dashed black lines 
are the band structure of phosphorene without strain for reference, (b) The schematic diagram of 
bilayer strained phosphorene-free phophorene heterostructure. The strain of the 1®* layer can be 
induced by the substrate, whereas the 2“'^ layer is unstrained because of the weak van der Waals 
interaction. 


phosphorene are shown in Supplementary Materials. Here, we choose the phosphorene 
with 2% compressed strain (along the armchair direction) as a donor. Its band structnre is 
shown in Fig. 6. For comparison, the band strncture of phosphorene without strain is also 
presented in the same figure. As can be seen, the position of CMB can be effectively tuned 
by the strain effect. 2% compressed armchair strain can shift the CBM down to the Fermi 
level around 0.11 eV because the compressed strain enhances the interaction of hybridized 
s — p orbitals of P atoms, which contribute to the CBM. The calculated PCE value of -2% 
strained-phosphorene/phosphorene can be around 20% as shown by the white solid circle in 
Fig. 5. Based on the theoretical studies of mechanical properties of monolayer phosphorene, 
the mechanical stability of phosphorene can be up to under 30% strain |20|, l60|. The 2% 
compressed strain of a monolayer phorphorene can be realized by the substrate effect in the 
experiment [see Fig. 6(b)]. Meanwhile, there is no strain on the second deposited monolayer 
phosphorene because of the weak van der Waals interaction between the two phosphorene 
layers. This can realize the strain/non-strain phosphorene heterostructure. 

Alternatively multilayer structures of phosphorene or TMDs may also be considered to 
increase the overall power conversion per unit area. Our calculated optical gap of phos¬ 
phorene is 1.6 eV which is at the edge of the infrared region. Therefore, heterostructures 
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with a phosphorene donor may be considered for multijunction solar cells absorbing pho¬ 
tons across the entire visible spectrum. Such cells incorporate several junctions that aim to 
absorb different portions of the solar spectrum so as to maximise total absorption [See Sup¬ 
plementary Materials] . Given the anisotropy of phosphorene, the stacking orientation in 
multilayer structures may be a significant factor 26|. 


IV. CONCLUSIONS 

In conclusion, through GW calculations, the band structures and optical spectra of mono- 
layer phosphorene have been calculated. The electronic gap of 2.15 eV and the optical gap 
of 1.6 eV are desirable for solar cell applications because of the strong exciton binding en¬ 
ergy. When paired with ZrS 2 or MoTe 2 the power conversion efficiency of the excitonic solar 
cells can be as high as 12%. There is further potential to improve the PCE of phospho¬ 
rene based solar cells substantially by tuning the materials, such as through the strain and 
multi-stacking effect, to achieve a better band alignment. 
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